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Abstract 


The  recently  developed  Laser  Raman,  Laser  Doppler  Velocimetry 
and  Coherent  Antistokes  Raman  Scattering  are  applied  to  the  diag- 
nostics of  flow  fields  and  flames.  The  concentration  of  species 
in  a cold  jet,  and  the  velocities  are  measured  and  compared  to 
measurements  using  standard  techniques.  The  Raman  and  LDV  tech- 
niques are  then  applied  to  diffusion  flames.  Experimental  results 
concerning  concentration  of  species,  temperature,  velocities  and 
turbulent  intensities  are  obtained  simultaneously.  The  latter  are 
obtained  from  the  LDV  data  as  well  as  the  concentration  data.  It 
is  seen  that  by  proper  processing  of  the  concentration  data  an 
indication  and  a measure  of  the  turbulent  intensity  may  be  obtained. 
It  is  further  shown  that  from  the  simultaneously  acquired  concentra- 
tion data  a parameter  of  major  importance  in  turbulence  modeling, 
the  "mixedness"  parameter  or  the  cross  correlation  function  may  be 
directly  measured. 
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TEMPERATURE,  CONCENTRATION 
VELOCITY  AND  TURBULENCE 
MEASUREMENTS  IN  JETS  AND  FLAMES 

S.  Lederman,  A.  Celentano  and  J.  Glaser 
1.  Introduction 

The  development  of  new  diagnostic  techniques  applicable  to 
fluid  dynamic  research,  has  been  an  ongoing  task  in  our  labora- 
tory for  many  years.  With  the  funding  by  Project  Squid  this 
task  was  directed  towards  optical  nonintrusive  techniques,  ap- 
plicable to  the  diagnostics,  of  flow  fields  propulsive  devices, 
and  combustion. 

As  a result  of  the  energy  crisis,  and  the  fact  that  almost 
all  of  the  energy  derived  from  fossil  fuels  is  obtained  as  a 
product  of  a combustion  process,  the  interest  in  the  detailed 
understanding  of  the  phenomena  involved  in  combustion,  has 
been  revitalized.  It  became  clear  that  a better  understanding 
of  the  processes  involved  is  vital  not  only  for  more  efficient 
designs  of  the  combustors,  with  the  associated  savings  in  fuels, 
but  from  the  environmental  point  of  view,  it  could  provide  an 
answer  to  the  problem  of  protecting  the  environment  by  a reduc- 
tion or  a possible  elimination  of  the  harmful  exhaust  emissions 
of  combustion  systems.  It  is  well  known  that  contrary  to  the 
highly  developed  technology  in  practical  design  of  combustion 
systems,  a scientific  understanding  of  the  details  of  the  com- 
bustion processes  is  still  in  its  infancy.  However,  recent  in- 
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vestigations  and  exchanges  of  views  in  the  scientific  community 
(Ref.  1,2, 3 ,4, 5)  have  contributed  greatly  to  a systematic  defini- 
tion of  the  problems  and  delineation  of  some  particular  aspects 
of  combustion  which  appear  to  be  of  major  importance  in  reaching 
a detailed  understanding  of  the  phenomena  involved.  It  has  been 
found  that  a phenomenon  exerting  a great  deal  of  influence  on 
the  combustion  processes  is  turbulence.  A knowledge  of  the 
turbulence  level,  instantaneous  temperature  and  concentration  of 
the  species  involved  may  be  very  useful  in  understanding  the  com- 
bustion process.  Modern  developments  in  nonintrusive  laser  diag- 
nostic techniques  of  flow  fields  could  be  of  significant  value  in 
this  context,  if  applicable  to  combustion  processes.  As  men- 
tioned previously,  nonintrusiveness  of  the  measurement  is  one  of 
the  major  requirements. 

It  has  been  generally  agreed  by  most  researchers  in  the  field 
that  the  first  aspect  of  nonintrusiveness  and  by  the  same  token 
noninterference,  can  be  best  met  by  optical  means.  Neglecting 
for  the  moment  some  of  the  difficulties  encountered  in  the  practi- 
cal implementation  of  some  of  the  optical  measurement  techniques, 
one  must  immediately  distinguish  between  the  internal  and  external 
flow  fields.  The  diagnostics  of  the  external  combustion  flow 
field,  having  almost  unlimited  optical  accessibility,  presents  no 
special  difficulties.  However,  the  diagnostics  of  the  flow  of 
the  internal  combustion  chamber  may  present  problems  not  only  of 
accessibility  but  also  of  potentially  disruptive  and  flow  dis- 
turbing features  introduced  by  the  necessary  modifications  of  the 
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combustion  chamber.  These  have  to  be  considered  and  minimized 
as  much  as  possible  in  each  individual  case. 

As  mentioned  above,  optical  methods  represent  the  best  hope 
at  present  of  achieving  the  best  measurement  in  combustion  with 
the  least  interference.  There  are  a number  of  optical  methods 
utilizable  for  diagnostic  purposes.  They  are  generally  associ- 
ated with  particle  and  molecular  scattering,  or  molecular 
absorption.  In  the  last  few  years  some  of  these  methods  have 
been  developed  to  the  point  where  they  are  now  incorporated  as 
part  of  commercially  available  instruments.  Most,  however,  are 
only  utilizable  at  present  in  research  laboratories.  The  major 
optical  methods  presently  utilized  for  diagnostic  purposes  are 
Mie  scattering  in  LDV,  Rayleigh,  in  density  and  temperature, 
fluorescence  in  concentration,  absorption  in  concentration, 
vibrational  Raman  in  concentration  and  temperature,  and  most  re- 
cently coherent  Antistokes  Raman  scattering  which  is  very 
promising  by  virtue  of  its  very  strong  signal  of  about  six  orders 
of  magnitude  higher  than  vibrational  Raman.  Most  of  the  above 
listed  diagnostic  methods  have  their  advantages  and  disadvantages. 
The  theoretical  background  of  most  of  the  above  scattering  phenomena 
are  different  and  require  generally  different  analytical  treatment. 
They  all  however  can  be  characterized  by  a common  convenient  param- 
eter known  as  the  equivalent  scattering  cross-section  in  units 
2 

cm  /Steradian.  This  parameter,  besides  being  a strong  function  of 
the  phenomenon  itself,  depends  on  the  character  of  the  scattering 
particle  and  the  frequency  of  the  illuminating  light,  among  others. 
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These  scattering  cross-sections  vary  over  an  extremely  large 

range.  Typically,  the  scattering  cross-section  for  the  Mie 
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scattering  is  in  the  range  of  10  cm  /Ster.  for  absorption  and 

-20 

fluorescence  of  the  order  of  10  , for  Rayleigh  of  the  order 

of  10  cm  /Ster.  and  Raman  of  the  order  of  10  cm  /Ster. 

CARS  may  have  a scattering  cross  section  of  up  to  6 orders  of 
magnitude  higher  than  the  regular  vibrational  Raman.  It  is 
obvious  that  a technique  with  the  larger  equivalent  scattering 
cross  section  is  to  be  preferred  over  the  ones  with  smaller 
cross  sections.  This,  however,  would  be  a great  simplification 
of  the  problem  at  hand  since  it  involves  only  one  of  the 
characteristic  parameters.  A number  of  criteria  have  to  be  met 
in  performing  a particular  measurement  in  a particular  situa- 
tion. Some  of  these  techniques  are  limited  in  sensitivity 
(small  scattering  cross  section)  but  are  capable  of  spatial  and 
time  resolution  (Raman) . others  may  be  more  sensitive  but  are 
incapable  of  spatial  or  time  resolution  (absorption) . Any  in- 
formation concerningspatdal  resolution  may  be  obtained  after  con 
siderable  accumulation  of  absorption  data  and  sophisticated  com 
putational  processing  of  the  same.  The  other  above-mentioned 
scattering  phenomena  except  the  Mie  scattering  are  generally 
superior  as  far  as  sensitivity  is  concerned  relative  to  Raman 
but  have  other  disadvantages  which  limit  their  usefulness.  A 
thorough  discussion  of  all  of  these  methods  is  beyond  the  scope 
of  this  report.  Since  in  the  last  few  years  in  the  course  of 
development  of  new  diagnostic  techniques  for  fluid  dynamic  re- 
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search  utilizing  lasers  our  efforts  were  directed  at  Laser  Raman 
Scattering  and  Laser  Doppler  anenometry,  this  paper  will  discuss 
some  of  the  results  of  our  efforts. 

Specifically,  our  attention  was  focussed  on  the  applicability 
of  Raman  scattering  towards  the  determination  of  specie  concentra- 
tion and  temperature,  among  others,  in  a mixing  axisymmetric  jet 
and  air-methane  flame.  This  particular  task  was  undertaken  with 
the  explicit  aim  of  demonstrating  the  feasibility,  proving  the 
applicability,  and  establishing  the  technology  of  remote,  in- 
stantaneous and  simultaneous  determination  of  specie  concentration 
and  temperature  of  the  individual  specie  in  a flow  field  without 
mechai  cal  probes.  In  this  report  after  a short  review  of  the 
spontaneous  Raman  effect,  the  LDV  technique  and  the  Coherent 
Antistokes  Scattering  technique,  an  attempt  is  being  made  to  de- 
termine temperature,  concentration,  velocity  in  coaxial  jets  and 
flames.  After  demonstrating  this  capability,  an  attempt  is  being 
made  to  determine  simultaneously  the  concentration  of  several 
species  of  interest  in  a flame,  their  individual  temperatures  as 
well  as  the  turbulence  intensity.  The  concentration  and  tempera- 
ture is  obtained  using  the  spontaneous  Raman  effect,  and  the 
turbulent  intensity  by  means  of  a L.D.V..  Using  the  concentration 
and  temperature  data,  an  attempt  will  be  made  to  extract  the 
flucuation  intensity  and  compare  it  with  the  intensity  as  obtained 
using  the  L.D.V.  techniques.  Furthermore,  some  data  concerning 
specie  concentration  in  a flame  are  discussed  using  the  coherent 
antistokes  Raman  scattering  method  which,  in  spite  of  some  dis- 
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advantages  as  compared  to  the  spontaneous  Raman  effect,  may  be  of 
major  importance  in  applications  concerning  combustion. 

2.  The  Raman  Effect 

The  Raman  Ef fect^' ^ is  the  phenomenon  of  light  scatter- 
ing from  a material  medium,  whereby  the  light  undergoes  a wave- 
length change  and  the  scattering  molecules  an  energy  change  in  the 
scattering  process.  The  Raman  scattered  light  has  no  phase  rela- 
tionship with  the  incident  radiation.  The  Raman  shifts  correspond 
to  energy  differences  between  discrete  stationary  states  of  the 
scattering  system.  Classically,  the  Raman  Effect  can  be  described 
as  the  modulation  of  the  scattered  light  by  the  internal  motions 
of  the  scattering  molecules.  In  this  kind  of  analogy,  the  Raman 
lines  would  correspond  to  the  side  bands,  and  the  Rayleigh  light 
to  the  carrier  frequency.  This,  of  course,  would  result  in  the 
Stokes  and  Anti-Stokes  lines  having  the  same  intensity,  which  is 
not  the  case.  Quantum  theoretically,  the  incident  photons  collide 
elastically  or  inelastically  with  the  molecules  to  give  Rayleigh 
and  Raman  lines,  respectively,  with  the  inelastic  process  much 
less  probable  than  the  elastic.  When  an  inelastic  collision  oc- 
curs with  the  incident  photon  furnishing  energy  to  the  molecule 
raising  it  to  a higher  energy  level,  the  scattered  photon  being 
of  lower  energy,  gives  rise  to  the  Stokes  line.  If  the  scattering 
molecule  gives  up  energy  to  the  impinging  photon  and  moves  to  a 
lower  energy  state,  the  scattered  photon  gives  rise  to  the  Anti- 
Stokes  line.  Since  the  Anti-Stokes  line  must  originate  in  mole- 


cules  of  higher  energy  level,  which  are  less  abundant  at  normal 
temperatures,  the  Anti-Stokes  lines  would  be  expected  to  be  much 
weaker  than  the  Stokes  lines.  The  process  of  light  scattering 
can  thus  be  visualized,  as  the  absorption  of  an  incident  photon 
of  energy  E by  a molecule  of  a given  initial  state,  raising  the 
molecule  to  a "virtual"  state,  from  which  it  immediately  returns 
to  a final  stationary  state  emitting  a photon  of  the  difference 
energy  between  the  two  states  and  incident  energy  E.  The  process 
is  illustrated  in  Figure  1. 

This  general  qualitative  behavior,  holds  for  the  vibrational 
as  well  as  rotational  transitions,  with  the  appropriate  selection 
rules,  which  tend  to  limit  what  appears  to  be  an  extremely  large 
number  of  possible  transitions  and  consequently  large  number  of 
Raman  lines. 

Since  for  the  purpose  of  this  work  the  vibrational  Raman 
scattering  is  of  direct  interest,  it  is  worthwhile  to  examine  the 
vibrational  Raman  response.  It  consists  essentially  of  three 
branches : 

a)  The  intense  Q branch  for  which  AJ  = 0 

b)  The  much  weaker  0 branch  for  which  AJ  = -2 

c)  The  much  weaker  S branch  for  which  AJ  = +2 
of  the  same  intensity  as  the  0 branch. 

It  can  be  shown  that  only  about  1%  of  the  total  vibrational 
intensity  resides  in  the  0 and  S branches  and  as  such  is  of  minor 
importance  as  far  as  the  present  application  is  concerned. 

In  quantitative  terms  the  scattered  intensity  may  be  written  as 


where  C is  a calibration  constant  of  the  system,  N is  the  number 
of  eatterers,  a is  the  equivalent  scattering  cross-section,  0. 
solid  angle,  and  > the  scattering  length. 

Equation  1 written  in  terms  of  the  scattered  signal  photons. 


becomes : 


n = 
s 


EoNCT-'nr\> 


and  in  terms  of  a voltage  signal 


V = 
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E No. -P . f).  t*  'll  .G.e.R 
o o q 
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where  Eq  is  the  incident  laser  energy,  and  -n  the  optical  and 
quantum  efficiencies  respectively,  G the  gain  of  the  photomulti- 
plier, e the  electron  charge,  R the  load  resistance,  E^  the  ener- 
gy of  the  scattered  photon,  and  t the  laser  pulse  duration.  In 
thermal  equilibrium,  the  ratio  of  the  Stokes  to  anti-Stokes  in- 
tensity can  provide  the  temperature  according  to  the  equation: 

- - 


It  is  clear  from  the  above  that  in  principle  it  is  possible  using 
Raman  scattering  to  obtain  instantaneously  and  simultaneously  the 
temperature  and  specie  concentration  in  a mixture  of  gases.  The 
former  because  the  Raman  transitions  take  place  in  a time  of  the 
order  of  fractions  of  pico-seconds,  the  latter,  because  both  the 
Stokes  and  anti-stokes  intensities  may  be  obtained  simultaneously. 

A major  drawback  of  the  Raman  diagnostic  technique,  is  the  ex- 
tremely small  equivalent  scattering  cross-section,  which  depend- 
ing on  the  incident  laser  frequency  and  specie  of  interest  may 
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vary  from  10  cm /sr  to  10  cm  /sr.  This  low  scattering  cross- 
section  forces  not  only  a minimum  limit  on  the  resolvability  of 
specie  concentration,  which  remains  relatively  high,  and  may  also 
limit  the  resolution  of  small  fluctuations  in  concentration  in 
a flow  field. 

A recent  new  development  in  Raman  Spectroscopy  may,  in  some 
cases,  improve  these  conditions.  This  new  development,  known  as 
CARS  ( Coherent-Anti-Stokes  Raman  Scattering) , has  been  shown  to 
have  an  equivalent  Raman  scattering  cross-section  of  up  to  6 
orders  of  magnitude  higher  than  the  spontaneous  Raman  Effect. 
Consequently,  specie  concentration  levels,  of  several  orders  of 
magnitude  lower  than  before  can  be  expected  to  be  resolvable. 


3.  The  CARS  Effect 

The  coherent  Anti-Stokes  Raman  Scattering  Effect  or 
11  12 

CARS  ' may  be  qualitatively  described  as  a process  by  which  a 
photon  v,  interacts  with  a tunable  photon,  v2  (Stokes  photon  of 
the  given  specie  of  interest)  through  the  third  order  non-linear 
susceptibility  to  generate  a polarization  component  of  the  Anti- 
Stokes  frequency  V3=2v^-v2.  This  is  diagrammatically  represented 
in  Fig.  2. 

Quantitatively  the  Anti-Stokes  scattered  power  can  be  shown  to 
be  represented  by: 
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where  P , P , P are  the  powers  of  the  Anti-Stokes,  incident 

3S  Li  S 

laser,  and  Stokes  radiation  respectively,  -tcQ^  is  the  coherence 
length  in  cm,  n is  the  index  of  refraction  at  the  Anti-Stokes 

9 S 

. . . . . 2 
frequency,  A is  the  interaction  cross-sectional  area  m cm 

X is  the  Anti-Stokes  wavelength  in  cm  and  x is  the  Raman 

9 S 

Susceptibility . 

The  coherence  length  -,coj1  defined  as  rr/ where  Ak=2k^-k2-k3 
may  be  written  as: 
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and  the  Raman  Susceptibility  y may  be  expressed  as 
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It  is  evident  from  Eq.  (5)  that  unlike  the  relation  of 
Eq.  (1),  the  specie  concentration  is  not  linearly  related  to  the 
scattered  radiation.  This  negative  feature  of  CARS  is  offset  by 
the  much  higher  equivalent  scattering  cross  section  (several 
orders  of  magnitude) , than  that  of  the  spontaneous  Raman  Effect. 

The  magnitude  of  the  equivalent  scattering  cross  section,  however, 
cannot  be  the  only  criterion  by  which  the  above  diagnostic  tech- 
niques may  be  evaluated.  Other  features  must  be  considered.  For 

example:  the  spontaneous  Raman  Effect  permits  the  measurement  of 

. 13 

many  species  which  may  be  present  in  a given  system,  simultaneous- 
ly using  a single  primary  laser.  This  is  not  possible  with  the  CARS 
diagnostic  method.  The  spontaneous  Raman  diagnostics  is  single 


ended 
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That  is,  the  transmitter  and  receiver  may  use  the  same 
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optics  or  may  be  located  in  proximity  to  each  other.  This  is 
not  possible  with  CARS  except  by  using  remotely  located  reflect- 
ing mirrors.  CARS,  on  the  other  hand,  due  to  its  coherent  high 
intensity  beam,  may  be  advantageous  in  systems  with  high  back- 
ground illumination,  fluorescence,  or  radiation  which  may,  in 
some  cases,  make  measurements  with  the  spontaneous  Raman  Effect 
impossible . 


4.  The  Laser  Doppler  Anemometer 
The  Laser  Doppler  velocitymeter , as  indicated  by  its  name, 
is  based  on  the  Doppler  principle.  Thus,  if  a small  volume  in 
a flow  field  is  illuminated  by  a laser,  the  frequency  of  the 
laser  light  scattered  by  moving  particles  in  the  volume  will 
appear  to  a receiving  stationary  photo-detector  as 
n V 

fn=f  ± t2-  . (e-e . ) (8) 

o 

where  fQ  and  are  the  frequency  and  wavelength  of  the  illum- 
inating laser  light,  V the  velocity  of  the  particle,  and  eg  and 
e^  the  scattered  and  incident  light  unit  vectors,  respectively. 

It  is  obvious  from  this  equation  that  in  principle  one  can 
measure  the  velocity  of  a particle  if  one  is  able  to  measure  the 
frequency  shift.  If  one  therefore  assumes  that  the  motion  of 
the  particle  whose  velocity  is  being  measured  is  equal  to  the 
motion  of  the  fluid  in  which  the  particle  is  immersed  the  motion 
of  the  fluid  can  be  measured.  This  assumption  can  be  made  only 
under  very  restrictive  conditions,  which  will  be  discussed  later. 


11 


As  with  the  Raman  scattering  technique,  the  Laser  Doppler 
velocitymeter  has  been  very  highly  developed  in  the  last  decade. 

A major  source  of  information  on  the  LDV  theory  and  operation  is 
References  15  and  16,  where  most  aspects  of  the  LDV  technology 
have  been  treated,  and  additional  references  can  be  found. 

In  contrast  to  the  Raman  scattering  technique,  which  is 
molecular  in  nature  and  is  essentially  omnidirectional  the  Laser 
Doppler  technique,  which  is  based  on  Mie  scattering,  is  highly 
directional.  In  other  words,  while  Raman  scattering  intensity 
under  certain  conditions  is  independent  of  the  angle  of  observa- 
tion, the  other  scattering  technique  is  highly  directional,  ex- 
hibiting a highly  pronounced  maximum  in  the  forward  direction. 

It  is  obvious  that  as  far  as  the  LDV  is  concerned  the  forward 
lobe  of  the  radiation  pattern  is  the  most  desirable  from  a variety 
of  points  of  view,  of  which  not  the  least  is  the  favorable  signal- 
to-noise  ratio. 

In  the  course  of  the  present  research  effort  in  our  laboratory 
concerning  laser  diagnostic  techniques  of  flow  fields,  the  prob- 
lem of  resolving  high  frequency  turbulence  spectra  was  looked  into. 
In  considering  the  LDV  for  this  problem  the  question  of  whether  or 
not  the  solid  particle  suspended  in  the  fluid  follows  the  stream- 
lines is  of  fundamental  importance.  How  far  do  the  results  of 
velocity,  as  well  as  that  of  frequency  spectra  of  tracer  particles, 
provide  the  requisite  information  regarding  the  turbulent  structure 
of  the  flow  field?  This  question  must  be  answered  through  careful 
analysis  and  consideration.  The  underlying  precept  permitting  LDV 


application  to  flow  field  diagnostics  is,  of  course,  the  assump- 
tion that  the  scattering  particles  are  moving  with  the  local 
fluid  velocity.  The  effect  of  particle  dynamics  on  the  perfor- 
mance of  an  LDV  have  been  considered  analytically  by  a number 
of  researchers  (Ref.  17-20)  in  the  field.  Yanta,  using  a gen- 
eralized drag  coefficient  examined  the  particle  dynamics  in  an 
expanding  supersonic  nozzle  and  noted  a lag  in  velocity  direct- 
ly proportional  to  the  diameter  of  the  scattering  particle. 

Tchen  and  Soo  examined  the  behavior  of  a particle  subjected  to 
homogenous  isotropic  turbulence.  Their  conclusions  were  that 
the  particle  diffusivity  was  the  same  as  the  Lagrangian  eddy 
diffusivity  of  turbulence.  However,  later  studies  by  Soo,  ex- 
perimental and  analytical,  based  on  the  "probability  of  en- 
counter" showed  that  these  two  dif fusivities  were  different  and 
that  the  seed  particles  do  not  generally  follow  fluid  particles. 
Even  for  very  small  particles  which  are  expected  to  follow  the 
flow,  the  increase  in  flow  Reynolds  number  decreases  the  particle 
diffusivity. 

In  Figure  3 the  ratio  of  the  particle  to  gas  velocity  as  a 
function  of  turbulent  frequency  with  the  particle  diameter  a_ 
as  a parameter  is  shown.  In  Figure  4 the  effect  of  2 particles 
simultaneously  present  in  the  scattering  volume  is  presented. 

It  is  evident  from  the  above  that  the  relationship  of  the 
frequency  spectra  from  an  LDV  to  the  turbulence  of  the  medium  is 
limited  to  the  lower  frequency  region,  that  limitation  being  a 


function  of  the  size  of  the  majority  of  the  scattering  particles. 
This  puts  also  a restriction  on  the  use  of  naturally  occurring 
particul.  es  because  their  sizes  are  generally  unknown.  The 
above  suggests  that  for  a more  meaningful  interpretation  of  the 
LDV  measurements,  a monitoring  of  the  sizes  of  the  particulates 
should  be  carried  out  and  incorporated  in  the  data  reduction 
process . 

Careful  choice,  however,  of  the  size  and  number  of  scatterers 
does  permit  one,  within  limits,  to  determine  the  velocity  of  the 
fluid  and  turbulent  intensity.  Since  the  major  problem  in  Laser 
Doppler  velocitymetry  is  the  acquisition,  processing,  and  hand- 
ling of  the  acquired  data,  not  the  principle  ifcself,  it  would 
serve  no  particular  purpose  to  go  into  the  different  optical 
arrangements  or  the  many  electronic  processing  methods  being 
utilized.  It  should,  however,  be  mentioned  that  the  dual-scatter 
type  which  can  be  operated  in  a forward,  as  well  as  a back  scatter 
mode  is  the  type  utilized  in  our  and  many  other  laboratories.  A 
schematic  diagram  of  the  dual  scatter  system  is  shown  in  Figure  5 

At  this  point  a few  remarks  are  in  order  as  far  as  the  dual 
scatter  system  is  concerned.  Due  to  the  symmetry  of  the  dual 
scatter  system,  the  output  current  of  the  photodetector  placed 
at  an  arbitrary  angle  with  respect  to  the  sample  volume  will  have 
an  AC  component  possessing  a frequency  proportional  to  the  parti- 
cles velocity  perpendicular  to  the  input  optics  axis  of  sym- 
metry (Up) . This  can  be  visualized  by  realizing  that  the  dual 
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beam  at  the  intersecting  point  generates  an  interference  fringe 
pattern.  Passage  of  a particle  through  successive  light  and 
dark  fringes  at  the  above  interference  pattern  will  result  in 
the  intensity  of  the  scattered  light  displaying  a sinusoidal 
variation  with  a frequency  equivalent  to  the  doppler  shift 
corresponding  to  U . Because  this  signal  is  now  independent  of 
the  direction  of  the  wavevector  for  the  scattered  radiation, 
light  may  be  collected  over  large  solid  angles.  The  subsequent 
increase  in  the  signal-to-noise  ratio  allows  this  mode  of  oper- 
ation to  be  used  in  situations  for  which  only  single  particles 
are  present  in  the  sample  volume. 

In  the  case  of  multiple  scatterers  present  at  the  same  time, 
the  sinusoidal  waves  mentioned  above  will  be  an  aggregate  of 
intefering  waves  as  indicated  in  Figure  4 and  will  result  in 
an  ambiguity  in  the  measured  velocity.  Therefore,  a sparingly 
seeded  flow  field  and  a properly  designed  electronic  processing 
system  is  recommended.  Furthermore,  each  situation  must  be  ex- 
amined for  optimal  results  and  least  error. 

Another  type  of  error,  due  to  the  statistical  biasing  of 
velocity  information,  has  been  examined  in  Refs.  21,  and  22. 

If  it  is  assumed  that  the  particulate  density  will  remain  con- 
stant throughout  the  flow,  the  rate  at  which  particles  pass 
through  the  sample  volume  will  be  a linear  function  of  velocity. 
This  will  cause  a biasing  of  the  measured  velocity  distribution 
towards  higher  velocities,  since  particles  traveling  at  these 
speeds  will  have  a larger  than  normal  probability  of  passing 


through  the  volume.  Tiederman,  Ref.  21,  has  considered  the  situa- 
tion in  which  the  true  velocity  distribution  is  gaussian  and  the 
u and  v components  are  totally  correlated.  The  results  indicated 
that  for  distributions  possessing  a standard  deviation  of  under 
15%  of  the  mean,  errors  in  the  measured  mean  velocity  were  under 
2%. 

At  this  point  it  should  be  noted  that  a simple  indication  of 
the  turbulent  intensity  may  be  obtained  using  an  LDV.  In  the 
case  of  a single  component  LDV  as  used  in  our  laboratories  de- 
fining the  turbulent  velocity  in  the  noraal  fashion  (u=u+u  ') , 
i.e.,  the  mean  value  in  addition  to  the  fluctuating  component, 
one  obtains. 


f = 2 si*1- 8/1  (u  + u')  (9) 

° Ao 

With  a large  number  of  velocity  samples  N taken  over  a peri- 
od of  time  large  with  respect  to  the  period  of  the  turbulent 
fluctuations,  it  is  possible  to  develop  a histogram  which  will 
represent  the  probability  distribution  of  the  velocity.  From 
this  it  is  possible  to  obtain  the  mean  velocity,  turbulence 
level  and  from  the  skewness  of  the  distribution  determine  if  the 
fluctuations  of  velocity  are  gaussian.  The  mean  velocity  will  be 


u = 


o^  1 1 
N2sinfi//2 


and  the  RMS  value  of  turbulence 


,2v*s  _ o 
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5.  Experimental  Apparatus 


The  experimental  apparatus  used  in  this  work  consisted  of 
4 basic  systems: 

a)  The  basic  Raman  System  used  to  obtain  temperature 
and  concentrations  of  a flow  field  or  flame. 

b)  The  L.D.V.  System  for  velocity  and  turbulent  in- 
tensity . 

c)  The  GARS  System  to  obtain  small  traces  of  unburned 
methane  in  an  air  methane  flame. 

d)  The  jet  and  combustion  system. 

A diagramatic  representation  of  the  complete  experimental  ap 
paratus  is  shown  in  Figures  6,  7 and  8 and  some  photographic 
views  in  Figures  9 and  10.  A full  description  of  the  experi- 
mental apparatus  is  given  in  Refs.  7,  13.  The  basic  components 
of  the  apparatus  are  indicated  in  the  schematic  diagrams. 

A common  component  in  the  above  experimental  arrangements  i 
the  data  processing  system,  consisting  of  a data  acquisition, 
storage  and  computing  facility.  This  system  allows  the  acquisi 
tion  of  a large  amount  of  experimental  data,  store  the  raw  data 
and  subsequently  retrieve,  analyze  and  present  in  a proper  man- 
ner. This  apparatus  permits  the  acquisition  of  data  concerning 
and  temperature  of  several  species  in  the  flow  field  as  well  as 
the  velocity  of  the  flow  simultaneously. 

6,  Experimental  Results 

Using  the  above  experimental  facility,  specie  concentration 
temperature  velocity  and  turbulent  intensity  of  a coaxial  jet 
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and  flame  have  been  obtained. 

Thus  Figure  11  presents  the  average  concentration  of  CO  2 
in  an  axisymmetric  jet,  as  obtained  using  a limited  number  of 
samples.  Figure  12  presents  a velocity  profile  on  the  same 
jet  using  several  measurement  techniques  and  comparing  the 
same.  Conceptually,  the  LDV  measurements  should  be  extremely 
accurate,  because  the  measurement  of  velocity  is  essentially 
reduced  to  the  measurement  of  frequency,  which,  as  is  well- 
known,  is  one  of  the  measurements  which  can  be  accomplished 
most  accurately.  However,  problems  in  seeding,  in  data  pro- 
cessing, etc.,  may  reduce  the  accuracy  of  the  velocity  being 
measured.  It  is  obvious  that  the  LDV  obtained  velocity  is  larger 
the  outer  portions  o^the  jet  than  those  obtained  using  alter- 
nate methods.  Part  of  this  difference  can  be  traced  to  the 
"statistical  biasing"  of  the  LDV  which  results  from  possible 
non-uniform  seeding  of  the  flow  field  in  particular  when  the 
flow  field  is  turbulent  in  nature,  part  from  the  nature  of  the 
electronic  processing  system  which  is  using  in  this  case  a high 
pass  filter,  thus  reducing  the  ability  to  measure  small  veloci- 
ties, and  partly  from  the  fact  that  the  scattering  particulates 
were  injected  into  the  jet  fluid  upstream  of  the  exit  orifice, 
and  no  effort  was  made  to  seed  the  surrounding  air.  Thus, 
towards  the  outer  edge  of  the  jet,  where  a large  portion  of  the 
fluid  has  been  entrained  from  the  surrounding  air,  the  particulate 
density  becomes  extremely  small.  (The  scatterers  may,  due  to 
inertia,  not  acquire  the  local  velocity)  and  the  acquisition  of 
data  becomes  tedious. 
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After  identification  of  the  possible  sources  of  error  in- 
volved in  the  measurements,  a series  of  tests  were  conducted 


on  a coaxial  jet  with  the  ratio  of  the  outer  to  inner  jet 

areas  (A  /A . ) = .26  and  U /U.  = .45.  The  results  are  shown 
o 1 o m 

in  Figure  13.  Figure  14  presents  the  turbulent  intensity 
corresponding  to  velocity  profiles  shown  in  Figure  13. 

In  Figures  15  and  16,  radial  distribution  of  velocity  pro- 
files for  A /A.  = 1.28  and  U /U.  = .5  and  U /U.  = .74  respective- 
o 1 o 1 o 1 

ly  are  shown.  A fairly  good  agreement  with  theory  and  other  ex- 
perimental data  is  readily  evident. 

Some  radial  distributions  of  velocity  in  a methane  air  flame 
are  shown  in  Figure  17  and  the  turbulent  intensity  in  Figure  18. 
Figure  19  and  20  present  examples  of  measurement  concentration 
and  temperature  profiles  in  an  air  methane  flame  respectively. 

The  above  velocity  concentration  and  temperature  data  have  not 
always  been  obtained  simultaneously.  However  the  recent  expan- 
sion of  our  data  acquisition  and  processing  systems,  permit  the 
routine  acquisition  of  such  data  simultaneously.  As  an  example 
of  this  new  capability  the  data  as  shown  in  Figures  21,  22,  23, 

24  and  25  have  been  obtained  simultaneously.  Thus  Figure  21 
presents  a velocity  profile  in  a flame  at  X/D  =5,2  with  the 
corresponding  computed  turbulent  intensity.  Figure  22  the  N2 
concentration  in  the  same  flame  at  the  same  X/D  and  the  corres- 
ponding concentration  fluctuation.  Figure  23  the  CO 2 concentra- 
tion and  concentration  fluctuation  and  Figure  24  and  25  the 
corresponding  N2  and  C02  temperatures  and  their  fluctuations. 


If  one  compares  the  N2  concentration  and  the  corresponding 
temperature  profiles,  somewhat  symmetrical  distributions  are 
evident.  This  is  not  the  case  for  C02 . Here  the  CC>2  con- 
centration is  higher  at  the  right  hand  side  of  the  distribu- 
tion profile,  while  the  corresponding  temperature  is  lower. 

This  is  perfectly  natural  since  cold  CC>2  was  introduced  into 
the  flame,  and  the  higher  C02  concentration  would  naturally 
correspond  to  a lower  temperature.  This  has  been  actually  ob- 
served visually  in  the  flame. 

As  a result  of  the  fact  that  the  data  were  obtained  simul- 
taneously, a parameter  of  importance  in  turbulence  modeling 
may  be  obtained.  This  parameter  referred  to  as  chemical 
’’mixedness"  is  shown  in  Figure  26  and  was  obtained  from  the 
concentration  data  of  Figures  21  and  23.  Finally  in  Figure  27, 
and  28  two  calibration  curves  for  Methane  and  Hydrogen  are 
given  as  obtained  using  the  CARS  apparatus  shown  in  Figure  7. 

An  attempt  to  obtain  the  traces  of  unburned  methane  in  a flame 
as  a function  of  axial  distance  from  the  exit  of  the  jet  using 
CARS  was  partially  successful.  However  the  ruby  laser  used  for 
that  purpose  has  broken  down  and  must  undergo  major  repairs. 

7.  Conclusions 

It  is  quite  clear  that  the  vibrational  Raman  scattering 
technique  and  the  Laser  Doppler  Velocimeter  are  capable  of  pro- 
viding nonintrusively  most  of  the  basic  information  regarding 
flow  fields.  The  limitations  of  those  techniques  appear  to  be 
more  in  the  practical  implementation  than  in  the  basic  concepts. 


Thus,  there  is  a lower  limit  of  resolution  capability  of  the 
Raman  scattering  technique  due  to  a combination  of  factors  such 
as  the  available  primary  laser  power,  the  number  of  scatterers 
in  the  sample  volume,  the  number  of  background  photons,  etc.. 

A very  convenient  parameter  to  assess  the  capability  of  a 
system  is  the  "feasibility  index".  Ref.  23.  This  index  was  de- 
fined as  X = NL-r  Qe  where  N is  the  number  density  of  the 
o 

scatterers  per  cm  , L is  the  length  of  the  sample  in  the  direc- 
tion of  the  laser  beam,  reference  cross  section,  and  0 and  e 
the  solid  angle  and  optical  efficiency,  respectively.  The  min. 
feasibility  index  for  a 1 j Ruby  laser  single  pulse  is  approxi- 
mately 10  Thus,  for  a situation  where  this  index  is  below 

10  15  a x joule  single  pulse  laser  would  not  provide  the  desired 
information.  An  increase  in  the  laser  energy  or  any  of  the 
other  factors  may  be  necessary.  There  is,  however,  a limit  on 
the  laser  energy  one  may  apply.  The  laser  energy  density  should 
be  below  the  breakdown  threshhold  which  for  Ruby  and  air  appears 
to  be  around  10  Watts/cm  . As  far  as  the  LDV  is  concerned,  it 
is  without  doubt  one  of  the  major  diagnostic  developments  for 
fluid  dynamic  research  attributable  to  lasers.  It  reduces 
velocity  measurements  to  the  measurement  of  frequency,  independent 
of  the  thermodynamic  state  of  the  medium,  capable  of  yielding  ac- 
curacies far  in  excess  of  any  other  method,  and  most  importantly, 
without  disturbing  probes  and  the  measurement  is  accomplished  re- 
motely. 


It  is  evident  from  the  above  that  the  optical  methods 
and  particularly  LDV  and  spontaneous  Raman  scattering  while  not 
free  from  difficulties  and  limitations,  are  essentially  the  two 
most  important  diagnostic  techniques  applicable  in  general  to 
flow  fields  as  well  as  to  external  and  internal  combustion. 

The  nonintrusive,  remote,  specific,  pointwise,  time  and  space 
resolution  capabilities,  are  the  major  assets  of  these  diag- 
nostic techniques.  The  low  scattering  cross  section  of  the  Raman 
diagnostic  techniqe,  which  is  the  source  of  major  difficulties, 
may  be  the  price  one  has  to  pay  for  the  other  advantages.  How- 
ever, improvement  in  the  collector  optics  and  data  acquisition 
equipment,  including  the  phtomultiplier  performance,  may  allevi- 
ate some  of  the  difficulties.  The  recent  developments  in  tunable 
dye  lasers  may  also  provide  some  possibilities  in  improving  the 
signal  to  noise  ratio  by  permitting  to  shift  into  a more  favorable 
frequency  band  in  terms  of  background  noise.  The  further  develop- 
ment of  the  (3^RS  method,  among  others  being  explored  at  the  moment, 
may  also  provide  a means  of  enhancing  the  optical  nonintrusive 
diagnostic  methods  in  flow  fields  and  combustion. 
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in  a cold  jet,  and  the  velocities  are  measured  and  compared  to 
measurements  using  standard  techniques.  The  Raman  and  LDV  tech- 
niques are  then  applied  to  diffusion  flames.  Experimental  results 

concerning  concentration  of  species,  temperature,  velocities  and  _ 
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20.  Abstract  (Continued) 

turbulent  intensities  are  obtained  simultaneously.  The  latter  are 
obtained  from  the  LDV  data  as  well  as  the  concentration  data.  It; 
is  seen  that  by  proper  processing  of  the  concentration  data  an 
indication  and  a measure  of  the  turbulent  intensity  may  be  obtained 
It  is  further  shown  that  from  the  simultaneously  acquired  concen- 
tration datp  a parameter  of  major  importance  in  turbulence  model- 
ing, the  <rmixedness **'' parameter  or  the  cross  correlation  function 
may  be  directly  measured. 
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